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1.0  Introduction 

Concrete  and  engineering  ceramics  exhibit  similarities  at  both  the  microscopic  and 
macroscopic  levels.  They  are  both  covalently  bonded  polycrystalline  materials  with  yield 
strengths  and  ductilities  that  are  strongly  dependent  on  pressure.  At  low  pressures,  they  fail  in 
brittle  fracture  at  relatively  small  loads;  and  at  high  pressures,  they  remain  elastic  to  much  larger 
loads,  and  incur  relatively  large  inelastic  strains  before  failure.  In  addition,  the  mechanical 
behavior  of  concrete  is  known  to  depend  upon  the  third  invariant  of  the  stress  tensor  [1-3]. 
Although  engineering  ceramics  are  more  homogeneous  and  less  porous  than  concrete,  the 
mechanical  and  micro  structural  similarities  of  these  two  classes  of  materials  raise  the  question 
whether  the  mechanical  behavior  of  ceramics  is  also  dependent  upon  the  third  invariant  of  the 
stress  tensor.  This  report  examines  the  possibility  that  the  yield  surface  of  silicon  carbide,  a 
ceramic  commonly  used  in  armor  applications,  is  dependent  upon  the  third  invariant. 

Traditional  techniques  for  measuring  yield  stresses  at  high  pressures,  such  as  tri-axial 
loading,  are  limited  by  the  very  high  strength  of  silicon  carbide  at  high  pressures.  As  a  result,  the 
best  means  of  evaluating  silicon-carbide  strength  characteristics  at  high  pressures  is  often  by 
high-velocity  impact.  In  this  report,  tests  of  high-velocity  impact  are  selected  from  the  literature 
for  investigating  the  influence  of  the  third  invariant.  The  tests  are  each  simulated  numerically 
with  several  values  of  the  tri-axial  yield  ratio,  a  material  parameter  defined  as  the  ratio  of  yield  in 
tri-axial  extension  to  yield  in  tri-axial  compression.  The  simulation  results  indicate  the 
sensitivity  of  each  test  to  silicon  carbide’s  dependence  on  the  third  invariant.  In  addition, 
comparison  of  the  test  and  simulation  results  suggests  the  actual  magnitude  of  silicon  carbide’s 
dependence  on  the  third  invariant. 

In  the  simulations,  the  JHB  ceramic  model  [4]  is  used  to  describe  the  mechanical  behavior 
of  silicon  carbide.  This  model  was  formulated  specifically  to  address  the  conditions  resulting 
from  high-velocity  impact:  large  strains  at  high  pressures  and  strain  rates.  Since  the  model  does 
not  include  the  influence  of  the  third  invariant,  it  has  been  augmented  for  the  simulations  in  this 
report  by  a  factor  dependant  upon  the  third  invariant.  Several  such  factors  have  been  proposed 
for  concrete  and  geomaterials,  e.g.  Mohr-Coulomb  [5],  Willam-Warnke  [6]  and  Gudehus  [7], 
each  of  which  depends  upon  the  tri-axial  yield  ratio.  Here,  the  factor  introduced  by  Gudehus  is 
incorporated  into  the  JHB  model.  In  addition,  several  of  the  original  JHB  yield  parameters  are 
re-calibrated  for  each  value  of  the  tri-axial  yield  ratio  in  order  to  maintain  agreement  with  the 
silicon-carbide  test  data.  This  process,  along  with  a  description  of  the  model  and  the  Gudehus 
factor,  is  described  in  the  next  section.  Simulation  and  test  results  are  presented  in  the 
subsequent  section,  followed  by  a  discussion  of  the  results.  Two  appendices  are  included 
documenting  the  influence  on  silicon  carbide  of  damage  softening  and  thermal  softening,  two 
additional  characteristics  for  which  the  JHB  model  has  been  modified. 
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2.0  The  JHB  Ceramic  Model 


2.1  Original  Form 

The  JHB  model  was  introduced  [4]  to  describe  the  behavior  of  ceramics  under  the 
conditions  of  ballistic  impact:  large  strains  at  high  strain  rates  and  pressures.  In  its  original  form, 
the  model  describes  the  yield  surface  as, 

a  =  (TTXC(p;D)(l  +  C\ns*)  (1) 

In  this  equation,  cr  is  the  equivalent  stress,  cr  =  ^3J2  ,  where  J2=s:s/ 2  is  the  second 

•  * 

invariant  of  the  stress  deviator,  s .  On  the  right  side,  p  is  the  pressure;  D  is  the  damage;  s  is 
the  dimensionless  strain  rate;  and  o^ip'fD)  is  the  isothermal  tri-axial  compressive  strength  at 

s  =1.0.  The  factor  (1  +  Cine  ) ,  including  material  parameter  C  ,  accounts  for  the  influence  of 
the  strain  rate.  Since  the  state  of  stress  is  only  represented  by  its  first  and  second  invariants 
(scaled  to  p  and  cr ,  respectively),  the  original  form  of  the  JHB  model  does  not  account  for 
effects  of  the  third  invariant. 

The  function  cr^ip^D)  takes  the  following  form. 


CT7JCC(p-,D)  = 


cr,  (P  +  T) 

( Pi+T ) 
Vi+iv max  -0-,.){l-eXP 


for  -T  <p<Pt 

-Vi(p-Pi) 

.( Pi  +  7’X°'max  "O’/). 


(2) 


for  p  >  pt 


This  equation  breaks  the  function  into  two  parts:  a  linear  function  for  —T  <  p  <  pt,  and  an 
exponential  function  for  p  >  pt .  The  two  functions  and  their  derivatives  are  continuous  at  their 
interface,  point  (p,,cr,.) .  The  model  parameters  are  pj ,  cr;,  the  asymptotic  limit  of  the  strength 
for  large  pressures,  crmax ,  and  the  tensile  (negative)  pressure  at  which  strength  vanishes,  T . 
These  parameters  are  shown  graphically  in  Fig.  1  for  the  yield  surface  of  undamaged  material, 
denoted  bycr7XC(p;D  =  0) . 

The  dependence  of  <jixc  (p:D)  on  damage  D  occurs  through  the  model  parameters.  In 
addition  to  the  set  of  parameters  depicted  in  Fig.  1  for  the  yield  surface  of  undamaged  material, 
another  set  defines  the  yield  surface  of  fully  damaged  material,  denoted  by  <jTXC(p\D  =  1). 

When  loading  reaches  the  undamaged  yield  surface,  the  material  incurs  plastic  strains  and 
damage  from  further  loading  each  timestep  n  . 

D  = 


z 


£Pf(P’K) 


(3) 


In  this  expression,  A s  (tn)  is  the  increment  in  plastic  strain  during  timestep  n,  and 
Spf  ( p,ln )  is  the  plastic  strain  to  failure  under  the  conditions  at  tn , 
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epf(P’tn)  =  Di 


r  p+t}n 

V  °"max  J 


(4) 


Here  Dx  and  N  are  model  parameters. 

Because  it  is  difficult  to  measure  damage,  there  is  some  uncertainty  about  the  strength  of 
partially  damaged  (0  <  D  <1 )  ceramics.  In  the  JHB  model,  the  strength  of  partially  damaged 
material  is  governed  by  the  undamaged  yield  surface.  When  damage  is  complete  (D  =1 ),  the 
strength  is  then  governed  by  the  damaged  yield  surface.  The  transition  between  these  two  yield 
surfaces  is  depicted  in  Fig.  1  by  the  dashed  line  indicating  instantaneous  softening  when  D  =  1 . 
An  alternate  approach  to  modeling  the  strength  of  partially  damaged  material  is  shown  by  the 
dashed  line  indicating  gradual  softening  when  D  >  0 . 


The  volumetric  response  of  the  ceramic  is  expressed  as, 

p  =Kj/r  +  K2ju2  +  K3p3  +Ap  (5) 


In  this  expression,  p  is  the  compressibility  (ratio  of  reduction  in  volume  to  initial  volume);  Apis 
the  increment  in  pressure  due  to  bulking;  and  Kl5  K2  and  K3  are  model  parameters. 


aTxc(P’D) 


Figure  1.  Dependence  of  tri -axial  compressive  yield  strength  on  pressure  in  the  JHB  ceramic 
model.  Undamaged  ceramic  is  governed  by  the  top  curve,  and  fully  damaged  ceramic  is  governed 
by  the  bottom  curve.  Model  parameters  pi?  ah  amax  and  T  are  shown  for  undamaged  ceramic. 


2.2  Incorporation  of  the  Third  Invariant 

Equation  1  expresses  the  yield  surface  in  terms  of  the  first  and  second  stress  invariants.  For 
the  incorporation  of  the  third  invariant,  principal- stress  space  provides  a  useful  context  because 
every  stress  state  can  be  uniquely  defined  by  its  three  invariants.  In  principal-stress  space,  the 
first  and  second  invariants  represent  the  line  cq  =  cr2  =  cr3  and  the  radial  distance  from  that  line, 

respectively.  The  portion  of  the  third  invariant  that  is  independent  of  the  first  and  second 
invariants  must  therefore  represent  the  angular  coordinate  of  a  cylindrical  coordinate  system  with 


4 
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axis  cr,  =  cr2  =  cr3  in  principal-stress  space.  The  Lode  angle  9  is  such  an  angular  coordinate, 
and  it  is  related  to  the  second  and  third  invariants  of  the  stress  deviators  ( J2  and  J3)  by  the 
relation  [8], 

j  ( 

sin  (3$)  =  —  —  (6) 

^  V^2y 

With  Lode  angle  9  as  the  independent  variable,  the  influence  of  the  third  stress  invariant  can  be 
introduced  into  the  JHB  model  by  appending  factor  Y(9:  T)  to  equation  1: 

cj  =  CTTXC(p;D)(l  +  an£*)r(9-x¥)  (7) 

Parameter  T  is  the  ratio  of  yield  strength  in  tri-axial  extension  to  yield  in  tri-axial 
compression,  and  the  factor  r(0;vP)  scales  the  tri-axial  compressive  strength  for  other  stress 
states.  The  form  of  T(#;VF)  introduced  by  Gudehus  [7]  is  used  in  this  study. 


1  +  sin  (30)  + 


'P 


(8) 


The  left  side  of  Fig.  2  shows  the  effect  of  r (9:  T)  on  the  yield  surface  in  principal-stress 
space  for  large  pressures  when  ¥  =  1.0,  T  =0.75  and  ¥  =  0.5.  The  view  is  of  a  plane  normal 
to  the  line  <j1  =  cr2  =  cr3 ,  in  which  Lode  angle  9  is  the  angular  coordinate,  and  the  curves  are  the 
intersections  of  the  yield  surfaces  with  the  plane.  The  curve  denoted  by  T  =  1.0  is  a  circle, 
indicating  that  T(0;vF  =  l.O)  =  1  for  all  9 ,  and  strength  is  independent  of  the  third  invariant. 
This  curve  represents  the  original  form  of  the  JHB  model,  as  Eqn.  7  reduces  to  Eqn.  1  when 
T  =1.  The  curves  denoted  by  T  =  0.75  and  T  =  0.5  demonstrate  successively  greater 
influence  of  the  third  invariant.  For  9  =  - 30°,  the  stress  state  is  tri-axial  compression 
(a,=cr2  >cr3),  T (-30°;^)  =  1.0,  and  the  curves  coincide  for  all  values  of  ¥  .  For  9  =  30°,  the 

stress  state  is  tri-axial  extension  ( al>cr2  =cr3),  r(30  ;T  )=  T ,  and  the  curves  differ 

maximally.  The  intermediate  case  of  9  =0°  represents  pure  shear. 
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Figure  2.  Effect  of  the  tri-axial  yield  ratio  'F  on  the  yield  surface, 
as  seen  in  a  plane  normal  to  the  line  a!  =  <t2  =  <r3  in  principal  stress  space 
for  high  pressures  (left)  and  low  pressures  (right).  0  is  the  Lode  angle. 


2.3  Parameters  for  Silicon  Carbide 

To  examine  the  influence  of  the  third  invariant  in  silicon  carbide,  simulations  of  plate- 
impact  and  ballistic  tests  are  compared  to  test  data  in  the  following  section.  The  simulations  are 
performed  with  'T  =1.0,  T  =0.75  and  T  =  0.5  to  identify  trends  in  the  influence  of  the  third 
invariant. 

The  yield  surface  of  undamaged  material  in  tri-axial  compression,  denoted 
by  a  rxc  ( p\D  =  0)  in  Fig.  1,  is  defined  by  parameters  pi ,  cr, ,  crmax  and  T  .  These  parameters  are 

unique  to  each  ceramic  and  they  are  generally  determined  by  best  agreement  with  compressive 
and  spall  plate-impact  and  Hopkinson-bar  test  data.  Since  each  of  these  tests  involves  high  strain 

rates,  the  term  (1  +  Cln£*)  in  Eqn.  7  must  be  factored  into  the  determination  of  the  parameters. 
In  similar  fashion,  the  term  170;  T)  must  be  factored  into  the  determination  when  T  <1.0,  i.e., 
the  influence  of  the  third  invariant  is  nonzero.  One  consequence  of  this  coupling  is  that  the 
parameters  that  best  agree  with  the  test  data  depend  in  general  upon  the  value  of  T  . 

Since  T  =  1 .0  for  all  values  of  T  when  the  state  of  stress  is  tri-axial  compression,  and  T 
varies  with  ¥  for  other  states  of  stress,  the  parameters  pi ,  crj ,  crmax  and  T  will  depend  upon 
T  only  when  they  are  derived  from  tests  in  which  the  state  of  stress  is  not  tri-axial  compression 
at  yield.  In  the  compressive  plate-impact  and  Hopkinson-bar  tests,  the  state  of  stress  at  yield  is 
indeed  tri-axial  compression.  Due  to  the  high  pressures  in  these  tests,  their  data  are  generally 
used  to  determine  the  exponential  part  of  Eqn.  2,  which  governs  when  p  >  pi .  However,  in  the 

spall  plate-impact  tests,  the  state  of  stress  at  yield  is  tri-axial  extension.  Since  the  pressure  is 
negative  (tensile)  in  these  tests,  their  data  are  used  to  determine  the  linear  part  of  Eqn.  2,  which 
governs  when  -T  <  p  <  pr  Therefore,  the  parameters  defining  the  linear  part  of  Eqn.  2  (  pt ,  cr; 
and  T )  must  be  redefined  for  each  value  of  T  to  ensure  that  the  model  agrees  with  the  spall 
tests.  But  the  parameters  defining  the  exponential  part  of  Eqn.  2  are  effectively  independent  of 
VF. 


6 


UNCLASSIFIED 


UNCLASSIFIED 


The  linear  parts  of  Eqn.  2  corresponding  to  each  of  the  three  values  of  ¥  are  depicted  by 
the  solid  lines  in  Fig.  3.  These  lines  represent  yield  in  tri-axial  compression  such  that  the 
corresponding  yield  in  tri-axial  extension  passes  through  the  spall  data.  These  lines  are  tangent 
to  the  same  exponential  curve  at  different  points.  Yield  in  tri-axial  extension  is  depicted  by  the 
dashed  and  dotted  lines  in  Fig.  3  for  ¥  =0.75  and  ¥  =  0.5,  respectively,  and  it  is  attained  by 
multiplying  the  yield  in  tri-axial  compression  by  r(30o;¥j  =  ¥ .  The  effect  of  redefining  the 
yield  surface  to  retain  agreement  with  the  spall  data  for  each  values  of  ¥  is  also  depicted  on  the 
right  side  of  Fig.  2,  where  the  yield  surfaces  intersecting  a  plane  normal  to  the  line  cr,  =  cr2  =  cr3 
can  be  seen  to  coincide  for  tri-axial  extension. 

In  the  simulations  reported  in  the  following  section,  the  JHB  model  parameters  used  for 
silicon  carbide  are  those  derived  by  Holmquist  and  Johnson  [9]  for  SiC-B,  except  for  the  values 
of  p,. ,  cr,  and  T  that  are  used  when  ¥  =  0.75  and  ¥  =  0.5 .  Those  values  are  shown  in  Table  1, 

along  with  the  values  of  pn  cr,  and  T  reported  by  Holmquist  and  Johnson  [9]  and  used  when 
¥  =  1.0. 


Table  1.  Values  of  silicon-carbide  yield-surface  parameters  that  depend  on  ¥ 


¥=1.0 

¥  =  0.75 

¥  =  0.5 

Pi  (GPa) 

1.50 

1.67 

2.01 

cr,(GPa) 

4.92 

5.28 

5.95 

T  (GPa) 

0.750 

0.876 

1.163 

Figure  3.  Redefinition  of  the  linear  part  of  the  JHB  tri-axial  compressive  yield  curve  for 
variations  in  ¥.  Redefinition  is  necessary  to  maintain  model  agreement  with  spall  data. 
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3.0  Computed  Results 

In  this  section,  the  results  from  simulations  of  high-velocity  impacts  with  silicon-carbide 
targets  are  reported  for  several  values  of  the  ratio  of  yield  in  tri-axial  compression  to  yield  in 
tiraxial  extension  (VF  =  1 .0 ,  lP  =  0.75  and  T  =0.5),  and  the  results  are  compared  to  test  data 
from  the  literature.  The  tests  include  tensile  and  compressive  plate  impacts,  long-rod  penetration 
of  several  thick  targets,  and  ballistic  impact  of  thin  targets  with  substrates.  Simulations  of  these 
tests  have  been  reported  by  Holmquist  and  Johnson  for  ¥  =  1.0  [9],  and  the  same  numerical 
algorithms  and  initial  finite-element  meshes  are  used  here.  Plate  impact  simulations  are 
performed  in  uniaxial  strain,  and  all  other  simulations  are  performed  in  two-dimensional 
axisymmetry.  Further  details  of  the  meshes  are  included  in  their  paper  [9]. 

3.1  Compressive  Plate  Impact 

Figure  4  shows  the  results  from  tests  and  simulations  of  two  compressive  plate  impacts:  at 
1566  m/s  on  the  top,  and  2259  m/s  on  the  bottom.  The  impacting  plate  is  silicon  carbide,  and  the 
target  plate  is  silicon  carbide  on  the  impacted  side  and  lithium  fluoride  on  the  back  side.  The 
plots  compare  the  computed  histories  of  silicon-carbide  velocity  at  its  interface  with  the  lithium 
fluoride  to  the  VISAR  test  results  of  Grady  and  Moody  [10].  Simulation  results  are  indicated  by 
a  line  of  long  dashes  for  'P  =1.0,  a  line  of  short  dashes  for  lP  =0.75 ,  and  a  dotted  line  for 
lP  =  0.5  ;  test  results  are  indicated  by  a  solid  line. 

Because  the  state  of  stress  in  the  silicon  carbide  is  tri-axial  compression  throughout  loading 
and  until  it  switches  to  tri-axial  extension  during  the  elastic  part  of  unloading,  the  three  computed 
results  are  identical  until  the  material  reaches  yield  in  tri-axial  extension.  Due  to  the  different 
values  of  lP  ,  the  computed  results  have  different  yield  points  in  tri-axial  extension.  This  is  seen 
by  divergence  of  the  three  computed  results  at  the  end  of  the  steeply  downward-sloping  part  of 
the  velocity  histories.  Yield  is  indicated  by  an  abrupt  transition  to  a  milder  downward  slope.  For 
both  impact  velocities,  the  results  from  'P  =  0.5  first  diverge  from  the  other  two  computed 
results  because  their  yield  in  tri-axial  extension  is  smallest.  The  results  from  'P  =  0.75  are  next 
to  make  the  transition,  as  they  represents  the  next  weakest  material  in  tri-axial  extension. 

In  comparing  the  computed  results  to  those  from  the  test  data,  the  transition  indicating  yield 
in  tri-axial  extension  is  best  replicated  at  both  impact  velocities  by  the  simulations  using 
lP  =  0.5  .  Further  examination  of  the  velocity  histories  after  this  transition  indicate  very  good 
agreement  between  the  simulations  using  ¥  =  0.5  and  the  test  data.  Agreement  with  test  data 
diminishes  for  the  simulations  using  T  =  0.75  and  'P  =  1.0 . 

In  these  simulations  of  compressive  plate-impact  tests,  pressure  remains  positive  as  the 
material  unloads  and  yields  in  tri-axial  extension.  As  a  result,  the  computations  compute  yield  in 
tri-axial  extension  by  factoring  the  exponential  part  of  the  yield  function  in  Eqn.  2.  Since  this 
part  is  the  same  for  each  of  the  computations  (independent  of  'P ),  yield  in  tri-axial  extension 
diminishes  with  diminishing  lP  ,  as  depicted  on  the  left  side  of  Fig.  2. 
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Figure  4.  Compressive  plate-impact  results  from  tests  [10]  and  simulations  for 
several  values  of  Y  and  impact  velocities  of  1566  m/s  (top)  and  2259  m/s  (bottom). 
Velocities  are  at  the  back  surface  of  the  silicon  carbide. 

3.2  Tensile-Spall  Plate  Impact 

Figure  5  shows  the  results  from  a  test  and  simulations  of  a  tensile-spall  plate  impact  at 
148  m/s.  Both  the  impacting  and  target  plates  are  silicon  carbide,  and  the  plot  compares  the 
computed  velocities  on  the  back  side  of  the  target  plate  to  the  VISAR  data  of  Dandekar  and 
Bartkowski  [11]. 

Unlike  the  simulations  of  compressive  plate-impact  tests,  the  pressure  in  tensile -spall  plate 
impact  tests  becomes  negative  as  the  material  unloads  and  yields  in  tri-axial  extension.  As  a 
result,  the  simulations  compute  yield  in  tri-axial  extension  by  factoring  the  linear  part  of  the  yield 
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function  in  Eqn.  2.  This  part  differs  for  each  value  of  ¥  in  order  to  maintain  the  same  yield 
stress  in  the  state  of  tri-axial  extension  at  spall,  as  shown  in  Fig.  3  and  the  right  side  of  Fig.  2. 
The  agreement  of  the  three  simulated  results  in  Fig.  5  is  therefore  confirmation  of  the  material 
parameters  of  Table  1,  as  they  were  derived  to  all  agree  with  the  spall  data. 


time  (ps) 


Figure  5.  Tensile-spall  plate-impact  results  from  the  test  [11]  and  simulations  for  several  values  of 

¥  at  an  impact  velocity  of  148  m/s.  Velocities  are  at  the  back  surface  of  the  silicon  carbide. 

3.3  Prestressed  Long-Rod  Penetration 

In  the  tests  of  Fundberg  et  al  [12],  silicon-carbide  targets  were  prestressed  by  confinement 
and  impacted  with  long  stationary  tungsten  rods.  Prestressing  was  achieved  by  heating  the  steel 
confining  tube  before  inserting  the  cylinder  of  ceramic.  Upon  cooling,  steel  plugs  were  screwed 
tightly  into  the  tube  ends.  Rings  were  then  welded  to  the  tube  ends  such  that  only  the  plugs  act 
as  cover  plates  along  the  target  axis.  During  the  tests,  x-rays  were  taken  at  several  times  to 
determine  penetration  depths. 

Figure  6  shows  a  comparison  of  silicon-carbide  penetration  versus  time  from  three  impact 
velocities.  The  test  results  are  indicated  by  squares  for  the  impact  velocity  of  1410  m/s, 
diamonds  for  1645  m/s,  and  triangles  for  2175  m/s.  The  simulations  were  performed  with  three 
values  of  ¥  for  each  impact  velocity.  They  are  represented  by  solid  lines  for  ¥  =1.0,  long- 
dashed  lines  for  ¥  =0.75,  and  short-dashed  lines  for  ¥  =  0.5.  For  the  impact  velocity  of 
1410  m/s,  the  rod  dwells  on  the  surface  of  the  silicon  carbide  in  the  test,  and  in  all  three  of  the 
simulations.  For  the  impact  velocity  of  2175  m/s,  the  rod  shows  negligible  dwell  and  penetrates 
at  a  nearly  constant  velocity  in  the  test  and  all  three  simulations.  For  the  intermediate  impact 
velocity  of  1645  m/s,  the  rod  dwells  before  penetrating.  Although  there  is  little  difference 
between  the  results  from  the  simulations  with  ¥=1.0  and  ¥  =  0.75,  the  simulation  with 
¥  =  0.5  transitions  from  dwell  to  penetration  about  3  ps  sooner  than  the  other  two.  The 
simulation  with  ¥  =1.0  was  used  by  Holmquist  and  Johnson  [9]  to  determine  the  damage 
parameters  D{  and  N  in  Eqn.  4.  So  for  these  simulations,  agreement  with  the  test  data  is  not  an 
indication  of  the  most  appropriate  value  of  ¥  . 
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Figure  6.  Results  of  tungsten  long-rod  impacts  with  thick  prestressed  silicon-carbide 
targets  from  tests  [12]  and  simulations  for  several  values  of  VF.  Test  data  are 
indicated  by  squares  (1410  m/s),  diamonds  (1645  m/s)  and  triangles  (2175  m/s). 


Figure  7.  Contours  of  sin(30)  at  10  ps  (left)  and  20  ps  (right)  from  the  simulation  of  a  long -rod 
impact  with  a  prestressed  silicon- carbide  target  at  1645  m/s  and  'F  =  0.5.  Values  of  sin(30)  near  tri- 
axial  compression  are  shaded  dark  gray;  those  near  tri-axial  extension  are  shaded  white. 

These  simulations  are  useful,  however,  in  demonstrating  the  influence  of  the  third  invariant 
on  the  very  sensitive  process  of  transitioning  from  dwell  to  penetration  of  a  ceramic.  The 
reduction  in  dwell  of  only  3  ps  between  'F  =  1.0  and  lF  =  0.5  indicates  little  influence.  Since 
the  pressure  is  positive  in  the  silicon  carbide  under  the  impact,  the  yield  surfaces  for  the  three 
values  of  *F  are  related  as  shown  on  the  left  side  of  Fig.  2,  such  that  yield  strengths  are  only 
similar  near  tri-axial  compression.  Therefore,  the  relative  lack  of  influence  of  'F  on  the  duration 
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of  dwell  indicates  that  the  silicon  carbide  under  the  impact  must  be  predominantly  in  tri-axial 
compression.  This  conclusion  is  verified  by  the  contour  plots  of  sin(30)  in  Fig.  7  from  the 
simulation  of  the  1645-m/s  impact  velocity  with  ¥  =  0.5  at  10  ps  and  20  ps.  In  these  plots, 
stress  states  near  tri-axial  compression  (sin (30)  <  -0.8 )  are  colored  dark  gray  and  they  are 
affected  very  little  by  'F  .  Stress  states  near  tri-axial  extension  ( sin(30  )>  0.8 )  are  colored  white 

and  they  are  strongly  affected  by  T* .  Intermediate  stress  states  are  colored  light  gray.  At  10  ps, 
the  rod  is  dwelling  and  all  of  the  silicon  carbide  under  the  impact  is  near  tri-axial  compression. 
At  20  ps,  the  rod  has  just  begun  penetration  and  the  loading  in  a  pocket  of  silicon  carbide  just 
below  the  impact  is  no  longer  near  tri-axial  compression. 


impact  velocity  (mA) 


Figure  8.  One-standard-deviation  error  bars  of  the  penetration  depths  from  simulations 
of  long -rod  impacts  with  confined  silicon -carbide  targets  for  several  values  of  T. 

Each  error  bar  represents  five  simulations,  and  two  of  them  are  shifted  horizontally 
at  each  impact  velocity  to  help  distinguish  them  visually.  The  error  bar  from  'F  =  1 
is  on  the  left  of  each  trio,  and  that  from  'F  =  0.5  is  on  the  right. 

3.4  Confined  Long-Rod  Penetration 

In  the  tests  of  Orphal  and  Franzen  [13],  deep  cylinders  of  silicon  carbide  were  confined  in 
titanium  tubes  with  aluminum  end  plates.  The  cylinders  were  then  impacted  with  long  tungsten 
rods  at  velocities  from  about  1500  m/s  to  just  under  5000  m/s.  In  this  section,  penetration  depths 
are  presented  from  simulations  of  the  tests  performed  at  impact  velocities  of  1300,  1500,  2000, 
3000  and  4000  m/s  for  the  three  values  of  lF  . 

Both  the  tests  and  the  simulations  exhibit  significantly  more  scatter  in  the  total  penetration 
than  in  the  primary  penetration,  apparently  due  to  the  highly  unstable  afterflow.  To  gauge  the 
influence  of  the  third  invariant,  comparisons  were  therefore  made  of  the  penetration  depths  at 
10  ps  after  impact,  a  time  when  primary  penetration  is  nearly  complete  and  no  afterflow  has 
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occurred  in  any  of  the  simulations.  As  a  result,  direct  comparison  to  the  test  data  is  foregone  in 
favor  of  less  scatter  in  the  results.  However,  this  does  not  elmimate  the  scatter,  as  long-rod 
penetration  of  confined  ceramics  has  been  shown  to  be  an  inherently  unsteady  process  [14]  in 
both  tests  and  computations.  Therefore,  five  simulations  were  performed  at  0.1-m/s  intervals  in 
each  of  the  five  impact  velocities  to  determine  the  sensitivity  of  the  penetration  depths  to  very 
small  perturbations.  For  example,  statistics  from  simulations  at  1999.8  m/s,  1999.9  m/s, 
2000.0  m/s,  2000.1  m/s  and  2000.2  m/s  represent  the  impact  velocity  of  2000  m/s. 

Figure  8  shows  the  penetration  depths  at  10  ps  normalized  by  the  original  rod  length,  L. 
The  five  simulations  performed  for  each  value  of  T  and  each  impact  velocity  are  represented  by 
a  one- standard-deviation  error  bar.  In  each  trio  of  error  bars,  two  are  shifted  horizontally  by  a 
small  increment  to  help  distinguish  them  visually,  with  the  error  bar  from  the  simulations  using 
¥  =  1.0  on  the  left,  that  from  T  =0.75  in  the  center,  and  that  from  T  =  0.5  on  the  right.  This 
plot  indicates  both  the  influence  of  the  third  invariant  and  the  sensitivity  to  small  perturbations  in 
the  impact  velocity.  For  each  of  the  impact  velocities,  the  mean  value  of  penetration  depth  (as 
indicated  by  the  midpoint  of  the  error  bar)  increases  with  decreasing  T* .  However,  for  impact 
velocities  of  2000  m/s,  3000  m/s  and  4000  m/s,  the  increases  are  small  compared  to  the 
magnitudes  of  the  penetration  depths. 

3.5  Unconfined  Long-Rod  Penetration 

In  the  tests  of  Behner  et  al  [15],  unconfined  silicon  carbide  was  impacted  onto  stationary 
gold  rods  at  velocities  of  3502  m/s,  3944  m/s  and  4462  m/s,  and  penetration  velocities  were 
determined  from  penetration  depths  recorded  by  x-rays.  For  this  study,  simulations  are 
performed  at  the  impact  velocity  of  4000  m/s  with  each  value  of  T* ,  and  the  penetration 
velocities  are  determined  analogously. 

Figure  9  shows  a  comparison  of  the  penetration  velocities  from  the  tests  and  simulations. 
The  effect  of  the  third  invariant  is  small,  and  each  of  the  simulations  agrees  well  with  the  test 
data.  Since  simulations  were  not  performed  with  perturbations  in  the  initial  conditions,  it  is  not 
possible  to  determine  their  sensitivities.  However,  there  is  no  clear  trend  in  the  influence  of  the 
third  invariant. 


impact  velocity  (m/s) 


Figure  9.  Penetration  velocities  of  long-rod  impacts  with  unconfined  silicon-carbide 
targets  from  tests  [15]  and  simulations  using  several  values  of  VF. 
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3.6  Thin  Target  with  Substrate 

In  the  tests  of  Wilkins  [16],  7.62-mm  pointed-nose  steel  bullets  were  shot  at  targets  of 
6.35-mm  silicon -carbide  tiles  over  6.35-mm  aluminum  plates,  and  a  ballistic  limit  of  660  m/s 
was  reported.  In  this  study,  simulations  of  these  impacts  are  performed  over  a  broad  range  of 
impact  velocities  at  4-m/s  intervals  for  each  value  of  . 


impact  velocity  (m/s) 


V  =  740  m/s,  t  =  10  ps,  ¥  =  0.5 


Figure  10.  Simulation  results  of  bullets  impacting  thin  targets  of  silicon  carbide  and  aluminum: 
bullet  residual  velocities  for  several  values  of  'F  (top)  and  contours  of  pressure  (bottom  left)  and 
sin(30)  (bottom  right)  from  the  simulation  with  VF  =  0.5  and  impact  velocity  of  740  m/s. 

The  top  of  Fig.  10  is  a  plot  of  the  residual  bullet  velocities  versus  impact  velocities  from 
each  of  the  simulations.  The  sensitivity  of  the  residual  velocities  to  variations  in  the  impact 
velocities  is  reflected  by  the  noise  in  the  data,  and  by  the  velocity  of  the  slowest  complete 
penetration  being  less  than  that  of  the  fastest  incomplete  penetration  for  both  ¥  =  1.0  and 
T  =  0.5  .  Despite  this  sensitivity,  the  data  show  a  clear  trend  toward  higher  ballistic  limits  with 
decreasing  T* ,  as  the  ballistic  limit  with  T  =  0.5  is  at  least  50  m/s  greater  than  that  with 
T  =  1.0.  This  trend  indicates  increasing  silicon-carbide  yield  strength  with  decreasing  'F  .  It  is 
therefore  the  opposite  trend  of  that  from  the  thick  targets  of  prestressed  and  confined  silicon 
carbide  of  preceding  sections.  Referring  to  Fig.  2,  increasing  silicon-carbide  yield  strengths  with 
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decreasing  *F  occur  only  in  the  linear  part  of  the  curve  (-T  <  p  <  pi ),  where  it  is  redefined  (via 
parameters  pn  cr,  and  T)  to  accommodate  the  spall  data.  The  contour  plots  on  the  bottom  of 
Fig.  10  verify  that  the  stress  states  in  the  thin-target  silicon  carbide  allow  it  to  access  these  higher 
yield  strengths  as  lF  decreases.  The  plots  are  from  the  simulation  with  lF  =  0.5  at  an  impact 
velocity  of  740  m/s  and  10  ps  after  impact.  For  this  value  of  *F  ,  yield  at  all  pressures  below 
Pi  =  2.0  GPa  is  governed  by  the  linear  part  of  the  yield  curve,  which  increases  for  decreasing 
*F  .  Pressures  below  pt  =  2.0 GPa  are  colored  light  gray  in  the  contour  plot  on  the  left.  Only  a 

small  region  of  the  silicon  carbide  just  beneath  the  impact  has  higher  pressure.  In  addition,  the 
plot  on  the  right  shows  contours  of  sin(3$) ,  and  it  indicates  that  most  of  the  silicon  carbide  is 
near  tri-axial  compression  or  intermediate  with  tri-axial  extension. 

The  silicon  carbide  tested  by  Wilkins  was  KT-SiC  [16],  not  the  hot-pressed  SiC-B  modeled 
in  these  simulations.  The  KT-SiC  was  of  lesser  density  and  lower  Hugoniot  elastic  limit,  and 
therefore  likely  of  lesser  ballistic  resistance.  As  a  result,  the  ballistic  limit  from  the  simulations 
should  be  greater  than  the  660  m/s  reported  by  Wilkins  [16].  In  the  plot  at  the  top  of  Fig.  10, 
only  the  ballistic  limit  from  the  simulations  using  lF  =  0.5  is  more  than  a  few  percent  greater 
than  660  m/s,  suggesting  that  'F  =  0.5  more  closely  represents  SiC-B  than  'F  =  1.0  or  'F  =  0.75  . 
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4.0  Discussion  and  Conclusions 

In  this  study,  the  Gudehus  [7]  form  of  dependence  on  the  third  invariant  is  incorporated  into 
the  yield  surface  of  the  JHB  ceramic  model.  The  model  parameters  for  silicon  carbide  are 
modified  to  retain  agreement  with  test  data  for  several  degrees  of  third-invariant  dependence. 
The  simulations  of  the  previous  sections  show  trends  in  the  influence  of  the  third  invariant  on  the 
behavior  of  silicon  carbide  under  high-velocity  impact. 

For  rods  impacting  thick  targets  of  silicon  carbide,  increasing  the  influence  of  the  third 
invariant  increases  penetration  depths  by  small  or  negligible  amounts.  In  prestressed  targets,  this 
effect  appears  as  a  small  decrease  in  the  duration  of  dwell  when  impact  velocities  are  in  the 
range  of  dwell  transitioning  to  penetration.  In  confined  but  non-prestressed  targets,  primary 
penetration  depths  increase  by  small  amounts  consistently  over  a  broad  range  of  impact 
velocities.  In  unconfined  targets,  the  third  invariant  shows  no  discemable  effect  on  the 
penetration  velocity. 

For  small-caliber  bullets  impacting  thin  targets  with  aluminum  substrates  (and  silicon- 
carbide  thickness  similar  to  bullet  caliber),  increasing  the  influence  of  the  third  invariant 
significantly  increases  the  ballistic  limit.  This  increase  in  performance  is  the  opposite  of  that  for 
thick  targets. 

In  addition  to  identifying  trends  in  the  influence  of  the  third  invariant,  it  is  also  possible  to 
infer  the  magnitude  of  the  influence  in  silicon  carbide.  This  can  be  done  by  comparing  the 
simulations  to  test  data.  In  the  compressive  plate-impact  tests,  reloading  to  yield  in  tension  is 
best  replicated  by  the  simulations  using  T  =  0.5  .  And  although  the  silicon  carbide  used  in  the 
ballistic  tests  of  thin  targets  is  not  the  SiC-B  simulated  in  this  study,  it  appears  the  simulated 
ballistic  limit  using  T  =  0.5  provides  best  fidelity  with  the  tests.  In  thick  targets,  the  influence 
of  the  third  invariant  is  too  small  relative  to  noise  in  the  tests  and  simulations  to  infer  the  most 
appropriate  influence.  Therefore,  it  appears  that  the  influence  of  the  third  invariant  is  significant 
in  silicon  carbide,  but  it  is  only  manifested  under  certain  loading  conditions  of  high-velocity 
impact,  such  as  those  in  thin  targets  and  compressive  plate  impacts. 
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Appendix  A:  Damage  Softening 

This  section  examines  the  effects  of  damage  softening  in  silicon  carbide  using  a  modified 
version  of  the  original  JHB  ceramic  model.  As  described  previously,  the  strength  of  partially 
damaged  material  (0  <  Z)  <  1)  is  governed  by  the  undamaged  yield  surface  in  the  original  JHB 
model.  When  damage  is  complete  (D  =  1 ),  the  strength  is  then  governed  by  the  damaged  yield 
surface.  As  depicted  in  Fig.  1,  the  transition  from  the  undamaged  to  the  damaged  yield  surface 
occurs  instantaneously. 

It  is  possible  that  the  strength  of  partially  damaged  material  is  not  governed  simply  by  the 
undamaged  yield  surface.  In  this  section,  two  alternate  descriptions  are  considered,  and 
simulations  are  performed  of  the  high-velocity  impact  tests  using  both  of  them.  In  the  first 
alternate  description,  the  yield  surface  is  the  undamaged  yield  surface  until  damage  is  halfway 
complete  (0  =  0.5);  subsequently,  it  is  interpolated  linearly  with  damage  between  the 
undamaged  surface  at  D  -  0.5  and  the  fully  damaged  surface  at  0  =  1.  In  the  second  alternate 
description,  the  yield  surface  is  interpolated  linearly  with  damage  between  the  undamaged 
surface  at  O  =  0  and  the  fully  damaged  surface  at  0  =  1.  This  description  is  depicted  in  Fig.  1 
by  the  dashed  line  labeled  “gradual  softening  when  O  >  0 .”  In  the  figures  showing  the 
simulation  results,  the  original  description  of  partially  damaged  yield  strength  is  denoted, 
“soft  =  0.0.”  The  first  alternate  description  is  denoted,  “soft  =0.5,”  and  the  second  alternate 
description  is  denoted,  “soft  =  1.0.”  With  this  nomenclature,  damage  softening  increases  as 
parameter  “soft”  increases. 

Figure  A-l  shows  the  simulation  and  test  results  of  the  compressive  plate  impacts.  For  both 
impact  velocities,  there  is  virtually  no  difference  between  the  velocity  histories  from  soft  =  0.0 
and  soft  =  0.5,  an  indication  that  negligible  damage  accumulates  beyond  0.5  in  the  silicon 
carbide.  However,  the  velocity  history  from  soft  =  1.0  differs  in  loading,  as  the  result  of  a  slower 
inelastic  wave  speed.  It  also  differs  upon  unloading  and  yield  in  tension,  as  the  pullback  signal  is 
smaller.  This  difference  brings  the  simulation  in  closer  agreement  with  the  test  data. 

Figure  A-2  shows  the  simulation  and  test  results  of  the  tensile- spall  plate  impacts.  Here 
there  are  insignificant  differences  between  the  three  degrees  of  damage  softening,  as  the  yield 
surface  is  not  reached  (and  damage  is  not  accumulated)  until  the  load  is  reversed  and  the  material 
fails  in  tension. 

Figure  A -3  shows  the  simulation  and  test  results  of  the  long-rod  impacts  with  prestressed 
targets.  Significant  differences  are  apparent  in  the  duration  of  dwell  for  the  1410-m/s  and  1645- 
m/s  impact  velocities.  In  particular,  at  1410  m/s  the  case  of  soft  =  0.0  fails  to  penetrate  after  at 
least  30  ps  of  dwell,  while  the  case  of  soft  =  1.0  dwells  for  only  about  5  ps.  At  2175  m/s,  none 
of  the  simulations  dwell  and  the  penetration  velocities  are  nearly  identical. 

Figure  A-4  shows  the  normalized  penetrations  at  10  ps  after  impact  from  the  simulations  of 
long-rod  impacts  with  confined  silicon-carbide  targets.  One  simulation  was  performed  at  each 
impact  velocity  and  each  method  of  damage  softening,  such  that  error  bars  indicating 
sensitivities  to  small  perturbations  are  not  incorporated.  The  data  show  increased  penetration 
depths  for  increased  damage  softening  at  impact  velocities  of  1300  m/s,  1500  m/s  and  2000  m/s, 
but  no  significant  differences  at  the  higher  velocities  of  3000  m/s  and  4000  m/s.  The  reduced 
strength  with  increased  damage  softening  is  apparently  only  significant  relative  to  inertial  effects 
at  the  low  impact  velocities. 
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Figure  A-l.  Compressive  plate-impact  results  from  tests  [10]  and  simulations  for 
several  degrees  of  damage  softening  and  impact  velocities  of  1566  m/s  (top)  and 
2259  m/s  (bottom).  Velocities  are  from  the  back  surface  of  the  silicon  carbide. 


Figure  A- 5  shows  the  simulated  and  test  penetration  velocities  due  to  long-rod  impacts  of 
unconfined  silicon-carbide  targets.  The  simulation  results  show  negligible  influence  of  damage 
softening  and  good  agreement  with  the  test  data. 

Figure  A-6  shows  the  residual  velocities  versus  impact  velocities  from  the  simulations  of 
small-caliber  bullets  impacting  thin  targets  of  silicon  carbide  with  aluminum  substrates.  The 
data  indicate  a  strong  influence  of  damage  softening  on  the  ballistic  limits,  with  the  ballistic  limit 
from  soft  =1.0  about  100  m/s  less  than  that  from  soft  =  0. 
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Figure  A-2.  Tensile-spall  plate-impact  results  from  the  test  [11]  and  simulations  for 
several  degrees  of  damage  softening  and  an  impact  velocity  of  148  m/s. 
Velocities  are  at  the  back  surface  of  the  silicon  carbide. 
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Figure  A-3.  Results  of  tungsten  long-rod  impacts  with  thick  prestressed  silicon- carbide  targets 
from  tests  [12]  and  simulations  for  several  degrees  of  damage  softening.  Test  data 
are  indicated  by  squares  (1410  m/s),  diamonds  (1645  m/s)  and  triangles  (2175  m/s). 


UNCLASSIFIED 


A-3 


UNCLASSIFIED 


1000  1500  2000  2500  3000  3500  4000  4500 


impact  velocity  (m/>) 

Figure  A-4.  Penetration  depths  from  simulations  of  long-rod  impacts 
with  confined  silicon-carbide  targets  for  several  degrees  of  damage  softening. 


3000  3500  4000  4500  5000 

impact  velocity  (nVs) 


Figure  A-5.  Penetration  velocities  of  long-rod  impacts  with  unconfined  silicon-carbide 
targets  from  tests  [15]  and  simulations  for  several  degrees  of  damage  softening. 
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Figure  A-6.  Residual  velocities  from  simulations  of  bullets  impacting  thin  targets 
of  silicon  carbide  and  aluminum  for  several  degrees  of  damage  softening. 


UNCLASSIFIED 


A-5 


UNCLASSIFIED 


A- 6 


UNCLASSIFIED 


UNCLASSIFIED 


Appendix  B:  Thermal  Softening 

This  section  examines  the  effects  of  thermal  softening  in  silicon  carbide  using  a  modified 
version  of  the  original  JHB  ceramic  model.  Thermal  softening  has  been  accounted  for  by 
appending  the  factor  onto  the  expression  of  the  yield  surface  in  Eqn.  1. 

<7  =  <JIXC(p',D)(l  +  Clns*)  (9) 

In  this  expression,  T  is  the  homologous  temperature,  T  =  (T  -Troom  )/(Tmelr  -  Troom ) ,  and 
m  is  a  model  parameter.  This  factor  reduces  the  strength  from  its  value  at  room  temperature, 
Troom ,  to  zero  at  the  melting  temperature,  Tmelt ,  and  its  form  has  been  adopted  from  the  Johnson- 

Cook  strength  model  for  metals  [17].  Indentation  tests  reported  by  Milman  et  al  [18]  of  SiC-B  at 
various  temperatures  have  been  used  to  find  the  model  parameters  for  the  simulations  in  this 
section:  m  =  2.6 and  Tmelt  =  1 1 40  AT .  However,  due  to  covalent  bonding,  the  yield  stress  of 
silicon  carbide  as  a  function  of  temperature  decreases  sharply  around  its  characteristic 
temperature  of  deformation,  above  which  crystal  lattice  resistance  to  dislocation  motion  is 
greatly  reduced.  In  addition,  the  function  shows  a  sharp  decrease  near  the  much  higher 
temperature  of  intensive  sublimation,  and  it  is  relatively  insensitive  at  other  temperatures.  The 
factor  does  not  allow  for  inflection  points,  and  the  thermal  softening  of  silicon  carbide 

can  therefore  only  be  accurately  modeled  to  about  600  degrees  C.  The  model  essentially 
neglects  all  strength  above  this  temperature. 

The  uncertainty  regarding  the  amount  of  internal- energy  growth  that  contributes  to  increased 
temperatures  has  led  to  the  inclusion  of  all  internal  energy  in  the  temperature  computations.  This 
assumption,  along  with  the  aforementioned  overestimate  of  thermal  softening  above  about 
600  C,  reduces  the  simulations  of  this  section  to  a  high  upper  bound  on  the  influence  of  thermal 
softening.  In  the  figures,  results  from  the  model  representing  the  upper  bound  on  thermal 
softening  are  denoted  by  m  =  2.6 ,  and  those  neglecting  thermal  softening  are  denoted  by  m  =  qo. 

Figure  B-l  shows  the  simulation  and  test  results  of  the  compressive  plate  impacts.  For  the 
impact  velocity  of  1566  m/s,  the  thermal  softening  results  in  a  significantly  reduced  inelastic 
wave  speed.  For  the  impact  velocity  of  2259  m/s,  thermal  softening  drastically  reduces  the 
inelastic  wave  speed,  and  also  reduces  the  magnitude  of  the  pullback  signal  upon  reloading  in 
tension. 

Figure  B-2  shows  the  simulation  and  test  results  of  the  tensile-spall  plate  impacts.  Since  the 
material  remains  elastic  until  it  reaches  the  yield  surface  after  reloading  in  tension,  there  are 
insufficient  inelastic  strains  to  activate  thermal  softening  before  failure. 

Figure  B-3  shows  the  simulation  and  test  results  of  the  long-rod  impacts  with  prestressed 
targets.  For  the  intermediate  impact  velocity  of  1645  m/s,  the  20  ps  of  dwell  are  completely 
eliminated  by  thermal  softening.  However,  the  duration  of  dwell  at  the  1410-m/s  velocity  is 
unaffected  by  the  thermal  softening,  and  the  penetration  velocity  at  the  2175-m/s  impact  velocity 
is  marginally  affected.  These  results  reflect  the  highly  unstable  nature  of  the  transition  from 
dwell  to  penetration. 
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Figure  B-l.  Compressive  plate- impact  results  from  tests  [10]  and  simulations 
with  and  without  thermal  softening  at  impact  velocities  of  1566  m/s  (top)  and 
2259  m/s  (bottom).  Velocities  are  from  the  back  surface  of  the  silicon  carbide. 

Figure  B-4  shows  the  normalized  penetrations  at  10  ps  after  impact  from  the  simulations  of 
long-rod  impacts  with  confined  silicon-carbide  targets.  One  simulation  was  performed  at  each 
impact  velocity  for  both  degrees  of  thermal  softening,  such  that  error  bars  indicating  sensitivities 
to  small  perturbations  are  not  incorporated.  As  for  damage  softening,  the  thermal  data  indicate  a 
clear  trend  toward  reduced  influence  of  softening  as  impact  velocities  increase.  For  the  impact 
velocity  of  1300  m/s,  penetration  occurs  only  with  thermal  softening;  and  for  the  impact  velocity 
of  1500  m/s,  penetration  increases  about  20  percent  with  thermal  softening.  For  the  higher 
impact  velocities,  the  influence  of  thermal  softening  becomes  negligible.  The  reduced  yield 
strengths  due  to  thermal  softening  are  apparently  only  significant  compared  to  inertial  effects  at 
the  lower  impact  velocities. 
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Figure  B-2.  Tensile-spall  plate-impact  results  from  the  test  [11]  and  simulations  with 
and  without  thermal  softening  at  an  impact  velocity  of  148  m/s.  Velocities  are 
at  the  back  surface  of  the  silicon  carbide. 

Figure  B-5  shows  the  simulated  and  test  penetration  velocities  due  to  long-rod  impacts  of 
unconfined  silicon-carbide  targets.  The  simulation  results  show  negligible  influence  of  thermal 
softening  and  good  agreement  with  the  test  data. 

Figure  B-6  shows  the  residual  velocities  versus  impact  velocities  from  the  simulations  of 
small-caliber  bullets  impacting  thin  targets  of  silicon  carbide  with  aluminum  substrates.  Unlike 
damage  softening,  thermal  softening  shows  no  influence  on  the  ballistic  limit.  The  difference  is 
perhaps  due  to  the  relatively  small  inelastic  strains  developed  before  failure  on  the  back  surface 
of  the  silicon  carbide.  Such  strains  are  necessary  to  invoke  thermal  softening,  while  damage 
softening  begins  immediately  upon  yield  for  soft  =  1.0. 
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Figure  B-3.  Results  of  tungsten  long- rod  impacts  with  thick  prestressed  silicon -carbide 
targets  from  tests  [12]  and  simulations  with  and  without  thermal  softening.  Test  data 
are  indicated  by  squares  (1410  m/s),  diamonds  (1645  m/s)  and  triangles  (2175  m/s). 
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Figure  B4.  Penetration  depths  from  simulations  of  long- rod  impacts  with 
confined  silicon- carbide  targets  with  and  without  thermal  softening. 
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Figure  B-5.  Penetration  velocities  of  long- rod  impacts  with  unconfined  silicon -carbide 
targets  from  tests  [15]  and  simulations  with  and  without  thermal  softening. 
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Figure  B-6.  Residual  velocities  from  simulations  of  bullets  impacting  thin  targets 
of  silicon  carbide  and  aluminum  with  and  without  thermal  softening. 
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